Solution of the Dirac Equation using the Lanczos Algorithm 
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Abstract 

Covergent eigensolutions of the Dirac Equation for a relativistic electron in an external Coulomb 
potential are obtained using the Lanczos Algorithm. A tri-diagonal matrix representation of the 
Dirac Hamiltonian operator is constructed iteratively and diagonalized after each iteration step to 
form a sequence of convergent eigenvalue solutions. Any spurious solutions which arise from the 
presence of continuum states can easily be identified. 
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INTRODUCTION 



In relativistic quantum mechanics, the calculation of the lowest lying bound states of a 
relativistic electron in a static Coulomb potential is a classic problem. In some respects 
the perturbative solution provided by Bethe[l| is both simple and elegant and a standard 
example in many textbooks. Early attempts to obtain the exact eigensolutions^, em- 
ployed finite difference methods to solve the radial equations and obtain the eigenenergies. 
Recently much effort has focused on the extension of basis set methodsj^, [s], 0, 0] which have 
been applied with success in non-relativistic quantum mechanics. Such matrix approxima- 
tions of the Dirac equation are usually obtained from a variational principle with the radial 
components approximated by a finite number of terms of a basis-set expansion. They often 
display pathological features such as the occurrence of unphysical spurious solutions 4, ^ 
and a failure to vary systematically with the basis-set size. These can for the most part be 
avoided if a resolvent operator rather than the Dirac Hamiltonian is used in the variational 
principle^] or by construction in explicit finite basis-set calculations [Qi]. Accurate numeri- 



cal solutions have also recently been obtained using mapped Fourier grid methods 
Spurious roots are still present in certain cases but can easily be identified jlo|. 



IG 



For the Coulomb potential with a coupling constant which is less than -^/|, the Dirac 



muii IS less iiifciii A / 

Hamiltonian is self-adjoint in the first Sobolev spacejl^. Hence, as we shall show, its 



eigenpairs can be obtained iteratively in an absolutely convergent manner by means of the 
Lanczos Algorithm 1^, Q provided, of course, that only the positive energy eigensolutions 
are considered. Although the Lanczos Algorithm is generally known for its application to 
conventional matrix eigenvalue problems, it can be used to solve eigenvalue problems of a 
self-adjoint operator which is bounded from below or above 1^. In the former case, the ma- 



trix to be diagonalized is iteratively transformed into a tri-diagonal matrix and diagonalized 
after each iteration step to form a sequence of convergent eigenvalue solutions. In the latter 
case, a tri-diagonal representation of the self-adjoint operator is constructed iteratively and 
diagonalized after each iteration step to form a sequence of convergent eigenvalue solutions. 



Any spurious solutions which arise 
tified in manner recently suggested 



Tom the presence of continuum states can easily be iden- 



151 ]. Furthermore, many of the aforementioned problems 



are avoided if the following iterative scheme proposed by Lanczos is employed. 
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THE LANCZOS METHOD 

Given the following eigenvalue problem 

^I^a) = ^aI^a) (1) 

where if is a self-adjoint operator (not necessarily bounded) in a separable Hilbert (or 
Sobolev) space which possesses a number of eigenvalues Ex (not necessarily bounded from 
above) which ascend from a minimal one Ei. Furthermore let |1) be a normalized start 
vector having the properties that: 

1. H^\l) exists for all non- negative integers k and 

2. satisfies the boundary conditions. 

The algorithm may then be stated as follows. Calculate successively the vectors 

\<l^i) ■■= |1) (2) 
1 " 

:= T, ^(i^ll) - J2 I0n')(0n'|i/"|1)), n = 1, 2, 3, . . . (3) 

where n is the iteration number. The Lanczos Method solves the eigenproblem by diagonal- 
izing successively the operators 

n 

Hn ■= ^ \(pm){4>m\H\4>m'){4>m'\, n = 1, 2, 3, . . . (4) 

m,m'=l 

The orthonormalized Lanczos vectors 10^) fulfill the following recursion relations 

|0„)=p„(^)|l) (5) 
where the Lanczos polynomials p„(x) are defined in the following way: 

pi{x):=l p2{x) := w^^{x - vi), (6) 

Pn+l{x) := W~^[{X - Vn)Pn{x) - Wn^lPn-l{x)] 72 = 2,3,4,..., (7) 

with 

Vn := {(f>n\H\(f)n) (8) 

Wn := {(l)n+l\H\(f)n) ■ (9) 



Performing the diagonalizations in equation H] reduces to finding the roots of the following 
characteristic polynomial 

det{Hn - X ■ i) := {-lYwi . . . (10) 

after each iteration step. 

The generated sequence of eigenpairs (exn , Wxn)) possess the following convergence 
properties 

exn ~^ Ex, A = 1,2,3,..., (11) 

n oo 



\exn) \Ex), A= 1,2,3,.... (12) 

n ^ oo 

or a self-adjoint operator the Lanczos algorithm can be implemented using Mathemat- 



ically 



IDENTIFYING SPURIOUS SOLUTIONS 



The exact bound states can be identified in the following manner 15|]. For those operators 
which possess a continuum as well as a point spectrum, the space spanned by the bound 
state eigenfunctions is by itself certainly not complete and a suitable start vector should 
be composed only of components in the subspace spanned by the bound state eigenvectors. 
Usually the start vector is chosen from a complete set of analytic functions which define a 
space J-'. This space is in most cases not necessarily of the same dimension as the subspace 
spanned by the exact eigenvectors. On the other hand, if the Lanczos algorithm is applied 
with this choice for the start vector, the eigenpairs obtained will correspond to those of 
the operator projected onto JF. A subset of these eigenstates must correspond to the exact 
eigenpairs of the unprojected Hamiltonian operator since the exact eigenstates can be ex- 
panded in terms of the complete set of states which span J-'. After each iteration, for each 
of the converging eigenpairs {exn, \e\n)), 

Axn = \eL - {exn\H^\exn)\ (13) 

(where n is the iteration number) is calculated and a determination is made as to whether 
A is converging toward zero or not. For the exact bound states of H, A must be identically 



zero while the other eigenstates states of the projected operator should converge to some 
non-zero positive value. Provided sufficient iterations are performed, it should be possible, in 
this manner to identify uniquely the approximate eigenpairs which ultimately will converge 
to the exact bound states. Since in the the Lanczos algorithm the lowest-lying eigenpairs 



(for an operator which possesses at most a lower bound) usually converge the quickest 171 ] 
one should be able to obtain good approximations to the lowest-lying eigenpairs. 

THE DIRAC EQUATIONS 

The coupled .adia, D.ac e,uatio„= fa the uppe. and lowe. co„.po„e„t. G(.) and F(.) 
can be written as|18l| 

_ hc^^ - ^^—-^F{r) = (E- V(r) - mc^)G{r), (14) 
dr r 

^JG{r)_ ^ (1 + i^)nc ^^^^ ^^^^ ^ rnc')Fir), (15) 

dr r 

where k is a quantum number analogous to the angular momentum quantum number / in 
non-relativistic quantum mechanics. Using f{r) = rF{r), g{r) = rG{r) and the Coulomb 
potential for charge Z, in relativistic units % = me = c = 1 the equations become 

-^ + -/(r) = (i?+--l)<?(r), (16) 
dr r r 

JiLl + _g^r) = {E+ — + 1 / r , 17 

dr r r 

where a is the fine structure constant. Equations f|T6|) and f|T7|) can written in matrix 
notation 

' <f-l)-ii-^)](^ir)\f,ir)\ 

Hd\^) = E\'^) (19) 

where 

l*)=f'"'*l (20) 

Equation (fT9!) defines the eigenvalue problem which will be solved iteratively using the 
Lanczos algorithm. 
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Iteration n 



FIG. 1: eAn VS. n for the bound state and the spurious solution 

NUMERICAL RESULTS 

To demonstrate that the the Lanczos algorithm may be used to solve the Dirac equations, 
a start vector was chosen using a slightly perturbed ground state solution derived by Bethe[l| 
with {k = —1) for an atom with z = 100. The start vector contained the following 
perturbed radial functions 

g^{r) = (r^)^e-""" (21) 

/i(r) = (22) 

za 

where s = \/l — (zaV and a perturbation factor 7 = .0998. Note this choice for the start 
vector eliminates the possibility of obtaining negative energy eigensolutionspl. 

Figure 1 shows the convergence for the ground state eigenvalue Con and a spurious solution 
ein- After 18 iterations, eo„ quickly converged to the exact ground state value of Eo =— 
0.683729. The unphysical spurious solution converged to a value of 3.37003. Using a larger 
perturbation factor of 7 = .95, eon converged to Eq to within two significant figures after 18 




FIG. 2: Aa n VS. n for the bound state and the spurious solution 

iterations. 

Figure 2 shows the convergence of Aa,„ for the ground state solution and the spurious 
solution. A\^n iteratively converges to zero for the generated bound state solution and 
converges to a non-zero value for the spurious solution. 



CONCLUSIONS 



The eigensolutions of the Dirac Hamiltonian have successfully been obtained using the 
the Lanczos algorithm. An eigensolution can be identified by testing the convergence of 
Ax^n for each generated eigenpair and by discarding generated solutions for which A does 
not converge to zero. The convergence rate, needless to say , depends on the choice of the 
initial start vector. 

Furthermore, it has been pointed out 17| that the Lanczos algorithm may be regarded as 
an application of the Rayleigh-Ritz method in which an orthonormalized set of Krylov vec- 
tors have been used to define the successive subspaces in which the operator is diagonalized. 
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This being the case, the convergence of A should also be a suitable means of sorting out the 
approximate exact eigenstates from the spurious eigenstates if the Rayleigh-Ritz method is 
employed. 
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